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Abstract 

Some recent experimental results obtained in collisions of heavy nuclei {y/s = 200 GeV) at 
BNL Relativistic Heavy-Ion Collider (RHIC) are discussed. The probes of dense matter created in 
heavy-ion collision by quarkonia, D and B mesons containing heavy charm and beauty quarks are 
considered. The centrality, rapidity and transverse momentum dependences of the nuclear modifi- 
cation factor and elliptic flow coefficient are presented and their possible theoretical interpretation 
is provided. 
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I. INTRODUCTION 



Lattice QCD (LQCD) calculations predict that at a critical temperature ^ 170 MeV, 
corresponding to an energy density Ec — 1 GeV/fm"^, nuclear matter undergoes a phase tran- 
sition to a deconfined state of quarks and gluons, called Quark-Gluon Plasma (QGP). At the 
modern collider facilities such as CERN supersynchrotron (SPS) (the nucleon-nucleon (NN) 
centre-of-mass energy for collisions of the heaviest ions is ^/s = 17.3 GeV), Relativistic Heavy 
Ion Collider (RHIC) at Brookhaven (^i = 200 GeV), and Large Hadron Collider (LHC) at 
CERN {y/s = 5.5TeV), whose heavy-ion program will start soon, heavy-ion collisions are 
used to attain the energy density, exceeding Ec- This makes the QCD phase transition po- 
tentially realizable within the reach of the laboratory experiments. The objective is then to 
identify and to assess suitable QGP signatures, allowing to study the properties of QGP. To 
that end, a variety of observables (probes) can be used jH-Q]. Further we will be mainly 
interested in heavy- flavour probes of QGP, i.e., utilizing particles having c- and b-quarks. 

A special role of heavy Q = (c, b) quarks as probes of the medium created in heavy-ion 
collision (HIC) resides on the fact that their masses {rric ~ 1.3 GeV, ^ 4.2 GeV) are sig- 
nificantly larger than the typically attained ambient temperatures or other nonperturbative 
scales, niQ ^ Tc,Aqcd = 0.2 GeV [5|. This has several implications: (i) The production 
of heavy quarks is essentially constrained to the early, primordial stages of HIC. Hence, 
heavy quarks can probe the properties of the dense matter produced early in the collision, 
(ii) Thermalization of heavy quarks is " delayed" relative to light quarks. One could expect 
that heavy quarks could "thermalize" to a certain extent, but not fully on a timescale of 
the lifetime of the QGP. Therefore, their spectra can be significantly modified, but still 
retain memory about their interaction history, and, hence, represent a valuable probe, (iii) 
RHIC, and especially LHC experiments allow to reach very low parton momentum fractions 
X, where gluon saturation effects become important. Heavy quarks are useful tools to study 
gluon saturation, since, due to their large masses, charm and bottom cross sections are 
calculable via perturbative QCD and their yield is sensitive to the initial gluon density. 

The heavy- flavor hadrons we will be interested in include: 1) open charm D = (cq) and 
open beauty B = {hq) mesons composed of a heavy quark Q = (c, b) and a light antiquark 
q = [u, d). These mesons could be sensitive to the energy density of the medium through the 
mechanism of in-medium energy loss; 2) hidden charm [charmonia= (cc)] and hidden beauty 
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[bottomonia=(66)] mesons (called collectively heavy quarkonia) being the bound states of 
the charm quark-antiquark, or bottom quark-antiquark pairs, respectively. Heavy quarkonia 
could be sensitive to the initial temperature of the system through the dissociation due to 
color screening of the color charge that will be discussed later. 

For detecting heavy flavor hadrons, different decay channels are used. At RHIC, in 
PHENIX and STAR experiments the measurement of the spectra of open heavy flavors 
is based on the measurement of the spectra of heavy flavor (HF) electrons and positrons 
[(e"*" + e~)/2] from the semileptonic decays like — K~e~^h'e, — K^e^Ve-, etc. These 
measurements are based on the fact that the decay kinematics of HF electrons/positrons 
largely conserves the spectral properties of the parent particles. Besides, the STAR ex- 
periment has the capability to directly reconstruct open heavy flavor mesons through the 
hadronic decay channels like D^{D^) — )■ i^^vr^, — )■ K^tx'^ti:'^ ^ etc. At low transverse mo- 
mentum, the STAR experiment also uses heavy flavor decay muons to provide open heavy 
flavor measurements. Quarkonia in both STAR and PHENIX experiments are detected 
through their dilepton decays QQ — )■ e^e~ (midrapidity) , QQ (forward rapidity). 



II. HEAVY FLAVOR PROBES OF QGP: QUARKONIA 

We begin the discussion of heavy flavor probes of QGP with heavy quarkonia. The 
question we would like to address is: What could happen with quarkonium yields in HIC if 
QGP is really formed? 

Let us note that if the deconfinement phase transition really takes place in a dense 
medium created in HIC, then a color charge in the QGP will be screened analogously to 
the Debye screening of an electric charge in the electromagnetic plasma. As a result of 
color Debye screening of the heavy quark interaction in QGP, the binding energy of a bound 
state decreases and one can expect that this would lead to the suppression of quarkonium 
yields in HIC. This idea was first suggested by H. Matsui and H. Satz {g] who predicted 
that color Debye screening will result in the suppression of J/ip meson (cc in the ^Si state, 
M = 3.097 GeV) yields. After that, the J/ip suppression was considered as one of the 
key probes for the QGP formation in heavy ion collisions. J/ip is especially promising 
because of the large production cross-section and dilepton decay channels which make it 
easily detectable. 
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However, soon it was realized that besides melting J/ip mesons in QGP due to the 
screening of the color charge, there are also a few competing mechanisms which could explain 
the suppression of J/il) production in heavy-ion collisions. These mechanisms are referred 
to as cold nuclear matter (CNM) effects. The first CNM effect is the absorption of J/il) by 
nuclear fragments from colliding nuclei. Let us consider, e.g., the proton- nucleus collision. 
Once produced in the hard primary parton processes, J/ip has to cross the length L of 
nuclear matter before exiting the nucleus, and, when traversing nuclear matter, it can be 
absorbed by forthcoming nucleons of a nucleus. The production cross section of J /ip in p-A 
collision can be parameterized as 



where app is the production cross section of J/ip in p-p collisions and o"^^ is the nuclear 
absorption cross section. From the global fit to the data on charmonium production in p-A 
collisions the value of a'^lf can be extracted, in particular, at SPS (NA50 experiment) it 
was obtained that a^^^ = 4.2 ± 0.5 mb 

The second CNM effect is related to shadowing of low momentum partons. This means the 
depletion of low momentum partons in nucleons bound in nuclei as compared to free nucleons. 
This effect can be accounted for in terms of the modification of the parton distribution 
functions in nucleon within the nucleus with respect to the parton distribution functions in 
a free nucleon: 



Here "i" denotes valence quarks, sea quarks, and gluons, x is the parton momentum fraction, 
is the momentum transfer squared. At high energies, J/ip^ are dominantly produced 
through the gluon fusion, and the J / ip yield is therefore sensitive to gluon shadowing. The 
underlying idea explaining the occurrence of gluon shadowing is that the gluon density 
strongly rises at small x to the point where gluon fusion, gg g, becomes significant. 
In the case of proton-nucleus and nucleus-nucleus collisions, where nuclei with large mass 
number A are involved, the nonlinear effects are enhanced by the larger density of gluons 
per unit transverse area of the colliding nuclei. A direct consequence of nuclear shadowing 
is the reduction of hard-scattering cross sections in the phase-space region characterized by 
small-x incoming partons. For gluons, e.g., shadowing becomes important at x < 5 x 10~^, 
and, hence, is relevant for the conditions of RHIC and LHC. Note however that the strength 
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(1) 





of the reduction is constrained by the current experimental data only for x > 10~^. 

Let us consider the J/ip production at RHIC experiments. The J/ip suppression can be 
characterized by a ratio called the nuclear modification factor 

d^Nff^JdpTdy 

obtained by normalizing the J/ip yield in A-B nucleus- nucleus collision by the J/ip yield in 
p-p collision at the same energy per nucleon pair times the average number of binary inelastic 
NN collisions. This ratio characterizes the impact of the medium on the particle spectrum. 
If heavy ion collision is a superposition of independent Ncoii inelastic NN collisions, then 
Rab = 1, whereas Rab < 1 {Rab > 1) corresponds to the case of the J/ip suppression 
(enhancement). As we discussed already, at first, it is necessary to clarify the role of CNM 
effects on J/t/.' production. At RHIC, CNM effects are studied in collisions of light deuteron 
and heavy gold nuclei, when the energy density reached in the collision is not enough for the 
formation of QGP. At RHIC energies, shadowing of partons is important and in the model 
calculations is implemented in two shadowing schemes for the nuclear parton distribution 



functions, the EKS model lOj and NDSG model The J/ip break-up cross sections 



obtained for two shadowing schemes 



and cTbroakup = 2.6^2:6 (NDSG) 



rom the best fit to data are o"brcakup = 2.8^2^ mb (EKS) 



12| . the term "break-up cross section" is 



3 (in Ref. 

used instead of the term "absorption cross section"). Although these values are consistent, 
within large uncertainties, with the corresponding value obtained at CERN SPS, a recent 



analysis [13[ shows that, in fact, the level of J/i/j CNM break-up significantly decreases with 
the collision energy. 

Let us now consider charmonium production in heavy ion collisions at RHIC. Figs. 1,2 
show the pT-iiitegrated J/i/j nuclear modification factor obtained in Au-Au collisions at 
RHIC/PHENIX experiment as a function of centrality, parametrized by the number of 



participating nucleons at mid- and forward rapidities, respectively |l2l . Il4| . The PHENIX 
data are shown by boxes. The Raa approaches unity for the peripheral collisions (small 
Npart) and goes down to approximately 0.2 at most central collisions (large Npart)- To see the 
level of the anomalous suppression beyond the cold nuclear matter effects it is necessary to 
extrapolate the CNM effects obtained in d-Au collisions to Au-Au collisions within the given 
shadowing scheme and corresponding J/ip break-up cross section. The results are shown by 
black and red curves with the corresponding error bands. It is seen that J/tp production 
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FIG. 1: (Color online) J/ip^s Raa for Au-Au collisions at midrapidity compared to a band of 
theoretical curves for the breakup values found to be consistent with the d-Au data. Both EKS 
and NDSG shadowing schemes are included. 




FIG. 2: (Color online) Same as in Fig. [T] but at forward rapidity. 

is significantly suppressed beyond CNM effects at forward rapidity and suppression is less 
pronounced at midrapidity in most central Au-Au collisions. 

Thus, one can conclude that assuming the conservative cold nuclear matter approaches 
some level of the anomalous suppression is indeed observed at RHIC which could characterize 
the produced medium as hot and deconfined. 
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FIG. 3: (Color online) J/^l' nuclear modification factor for the most energetic SPS (Pb-Pb) and 
RHIC (Au-Au) collisions, as a function of the number of participants Npart- 

However, not all is still clear. Measurements of the J/ip suppression by PHENIX col- 
laboration at RHIC lead to some surprising features. Fig. 3 shows compiled data for the 
nuclear modification factor obtained in CERN SPS and RHIC PHENIX experiments. There 
are two surprising results in these measurements. First, the mid rapidity suppression in 
PHENIX (the red boxes) is lower than the forward rapidity suppression (blue boxes) despite 
the experimental evidence that the energy density is higher at midrapidity than at forward 
rapidity, and, hence, one could expect that at midrapidity J /ip will be more suppressed due 
to higher density of color charges. Secondly, the nuclear modification factor Raa at midra- 
pidity in PHENIX (red boxes) and SPS (black crosses) are in agreement within error bars, 
a surprising result considering that the energy density reached at RHIC is larger than the 
one reached at SPS. This indicates that at RHIC energies additional mechanisms countering 
the suppression, could be operative. 

Let us consider possible explanations of the above features. 1. Regeneration of J/ip in 
the hot partonic phase from initially uncorrelated c and c quarks (quark coalescence model). 
If to compare the J/ip suppression pattern at RHIC and SPS, J/ip could be indeed more 
suppressed at RHIC than at SPS, but then regenerated during (or at the boundary of) the 
hot partonic phase from initially uncorrelated c and c quarks. If to compare the results at 
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FIG. 4: (Color online) J/tp^s V2 at midrapidity for [20, 60]% in centrality, as a function of pT, using 
42% of 2007 Au + Au statistics, with some theoretical predictions. 

RHIC (midrapidity vs. forward rapidity), at midrapidity, due to the higher energy density, 
there are more c and c quarks to regenerate than at forward rapidity that could explain the 
stronger suppression at forward rapidity. Note that the total number of initial cc pairs is 
larger than 10 in the most central Au-Au collisions. Certainly, if regeneration is important at 
the RHIC conditions, it will be even more important at the LHC conditions where more than 
100 cc pairs is expected to be produced in the central Pb-Pb collisions. 2. J/tp production 
could be more suppressed at forward rapidity due to the nuclear shadowing effects, which 
could be more pronounced away from midrapidity [isl ]. 

One of possible experiments aimed to verify the quark coalescence model is the measure- 
ment of the J/ elliptic flow |l^ . The idea is that if charmonia were produced by coalescence 
of charm quarks, they should inherit somehow their flow, which is known to be quite large 
from the open heavy flavor measurements (see the next Section), resulting in a higher V2 
than in the case of the direct production of J/ ip in hard collisions. The PHENIX experiment 
reported a first tentative measurement of J/tp^s V2 in Au-Au collisions As is seen from 
Fig. 4, these proof-of-principle measurements at the current level of precision do not allow 
one to distinguish between models assuming various level of regeneration (and thus elliptic 
flow) and much larger statistics is probably needed to differentiate between different models. 

To show the complexity of the problem, let us consider some theoretical models for 
the charmonium production, which quite satisfactory describe the RHIC data but whose 
predictions for LHC are drastically different. First, in the statistical hadronization model 
(SHM) 18[, it is assumed that: 1. All heavy quarks (charm and bottom) are produced in 
primary hard collisions and their total number stays constant until hadronization. 2. Heavy 
quarks reach thermal equilibrium in the QGP before the chemical freeze-out (hadronization). 
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FIG. 5: (Color online) Rapidity dependence of R^l^^ for two centrality classes in the statistical 
hadronization model. The data from the PHENIX experiment (symbols with errors) are compared 
to calculations (lines, see text). 

3. All quarkonia are produced (nonperturbatively) through the statistical coalescence of 
heavy quarks at hadronization. Multiplicities of various hadrons are calculated with the 
grand canonical ensemble. The generation of J proceeds effectively if c, c quarks are free 
to travel over large distances implying deconfinement. 

Fig. 5 shows the rapidity dependence of the nuclear modification factor, obtained in this 
model and the comparison with the rapidity dependence at PHENIX for two centrality bins. 
Two theoretical curves correspond to two fitting procedures, with one and two Gaussians of 
the J / ip data in pp collisions. In both cases, calculations reproduce rather well (considering 
the systematic errors) the Raa data. The model describes the larger suppression away from 
midrapidity. The maximum of Raa at midrapidity in this model is due to the enhanced 
generation of charmonium around midrapidity, determined by the rapidity dependence of 
the charm production cross section. The centrality dependence of Raa at ?/ = is shown 
in Fig. 6. The model reproduces quite well the decreasing trend with centrality seen in the 
RHIC data. Fig. 6 also shows the prediction of the model for the LHC. At much higher 
LHC energies, the charm production cross section is expected to be larger by about an order 
of magnitude. As a result, a totally opposite trend as a function of centrality is predicted, 
with Raa exceeding unity for central collisions. 



Let us consider the comovers interaction model (CIM) 19[]. This model does not assume 
the deconfinement phase transition. Anomalous suppression of J/ip (beyond CNM effects) 
is the result of the final state interaction of the cc pair with the dense medium produced 
in the collision (comovers interaction). The model consistently treats the initial and final 
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FIG. 6: (Color online) Centrality dependence of the R'K^ at midrapidity, according to the statistical 
hadronization model. 

state effects. The initial state effects include: 1) nuclear absorption of the pre-resonant cc 
pairs by nucleons of the colliding nuclei, 2) consistent treatment of nuclear shadowing for 
hard production of charmonium. The final state effects include absorption of the cc pairs 
by the dense medium created in the collision (interaction with comoving partons or hadrons 
produced in the collision). The model does not assume thermodynamic equilibrium and, 
thus, does not use thermodynamic concepts. The density of charmonium is governed by the 
differential rate equation 



dn 



T- 



dr 



-(^co[nco{b, s, y)nj/^{b, s, y)- 



(4) 



supposing a pure longitudinal expansion of the system and boost invariance. In Eq. (j4]), Uco 
is the density of comovers, which is found in the dual parton model 20| together with the 
proper shadowing correction, aco is the cross section of J/ ip dissociation due to interactions 
with comovers, taken such as to reproduce the low energy SPS experimental data (with 
(Jco = 0.65 mb). The first term on the right describes dissociation of charmonium due to 
interaction with comovers. The second term describes the recombination of charmonium and 
is proportional to the product of densities of charm quarks and antiquarks. The important 
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FIG. 7: Results for J /tjj suppression in Au-Au collisions at RHIC at midrapidity in the comovers 
interaction model. The solid curve is the final result. The dash-dotted one is the result without 
recombination (C = 0). 

feature of the CIM is that recombination of c-c quarks proceeds only locally, when the 
densities of quarks and antiquarks are taken at the same transverse coordinate s. This is 
different from the recombination in the SHM, where recombining quarks can be separated by 
large distance that implies deconfinement. The effective recombination cross section in the 
CIM is equal to the dissociation cross section due to the detailed balance. The results for the 
centrality dependence of the Raa are shown in Fig. 7. At midrapidity, the experimental data 
are well reproduced by full theoretical calculations (solid curve) taking into account nuclear 
shadowing, dissociation by comovers and recombination from charm quark and antiquark 
pairs. The results at forward rapidity, presented in Fig. [HI also well agree with the data, in 
particular, the J / ip suppression at forward rapidity is somewhat larger that the suppression 
at midrapidity. 

Fig. 9 shows the predictions of the model for LHC The parameter C encodes the recom- 
bination from c-c pairs and vanishes in the absence of recombination. Although the density 
of charm grows substantially from RHIC to LHC, the combined effect of initial-state shad- 
owing, absorption and comovers dissociation overcomes the effect of parton recombination. 
This is in sharp contrast with the predictions of the statistical hadronization model where 
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FIG. 8: Same as in Fig. 7, but at forward rapidity. The dashed Une is the total initial-state effect. 
The dotted line is the result of shadowing. In Fig. 7, the last two lines coincide. 

a strong enhancement of the J/ip yield with increasing centrahty was predicted. 

Thus, the J/ip suppression is an important characteristic to search for QGP. But if J/ if) 
anomalous suppression (beyond CNM effects) was observed in HIC, there are a few com- 
peting mechanisms to explain that: 1. Charmonium is dissociated due to the genuine color 
screening in the deconfined medium. 2. Charmonium is dissociated through interactions 
with comoving partons or hadrons in the medium formed in HIC. How to differ these mech- 
anisms? In comovers interaction model, the anomalous suppression sets in smoothly from 
peripheral to central collisions rather than in a sudden way when the corresponding disso- 
ciation temperature in the deconfined medium is reached. Even at SPS, where the role of 
quark recombination is of minor importance, current experimental errors still do not allow 
to disentangle these two mechanisms. However, even if the color screening mechanism is 
dominating, it is unclear what is really melted, directly produced J / ^/'s or originating from 
the feed-down of less bound charmonium states, Xc (— ^ J li^ + i^' J/i^ + which 
have lower dissociation temperatures. At RHIC and LHC conditions, the feed-down from B- 
meson decays becomes also important. Recombination enhances J ftp production and much 
complicates the picture but its effect may be different depending on whether the deconfine- 
ment phase transition happened or not. It is even possible that after all J/tps were melted 
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FIG. 9: Results for J/ij; suppression in Pb+Pb at LHC {^/s = 5.5 TeV) at midrapidity for different 
values of the parameter C in the comovers interaction model. The upper line is the suppression 
due to initial-state effects. 

in QGP they can be statistically regenerated at hadronization. True operative mechanisms 
of J/ production can be established only after studying all important dependences (from 
central! ty, rapidity, collision energy y/s,...) of all relevant observables with sufficient accu- 
racy. Further we will consider also the J/ip production in RHIC experiments as a function 
of the transverse momentum, but, for the sake of comparison, this will be more illustrative 
to do after presenting the respective dependence for open heavy flavor mesons. 

So far we considered the suppression (or enhancement) patterns for the J/ip produc- 
tion. Now let us briefly discuss the perspective for bottomonia. One could expect that 
bottomonium production might be easier to understand than charmonium production due 
to the following reasons. 1. Since less than one bb pair is produced in one central Au-Au 
collision, the regeneration is negligible at the conditions of RHIC. Besides, only about 5 
bb pairs are expected to be produced in a single central Pb-Pb collision at LHC. Hence, 
regeneration should play much less role in the beauty sector than in the charm sector. 2. 
Having higher masses, bottomonia originate from higher momentum partons and will less 
suffer from shadowing effects. 3. The absorption cross section for T is by 40 — 50 % smaller 
than the corresponding cross section for J/ip and 
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These features should ease the separation of the anomalous suppression in the T's family. 



III. HEAVY FLAVOR PROBES OF QGP: OPEN HEAVY FLAVOR MESONS 

What qualitative effects could one expect to obtain when probing the dense matter by 
heavy quarks (charm or bottom)? As well known from electrodynamics, the bremsstrahlung 
off an accelerated heavy quark Q is suppressed by the large power of its mass ~ {niq/mQ)^ as 
compared to light quarks. Therefore, gluon radiation off heavy quarks (i.e., radiative energy 
loss) is much suppressed relative to light quarks. Consequently, one could expect a decrease 
of high pt suppression and of the elliptic flow coefficient V2 from light to charm to bottom 
quarks. Or, that the energy loss and coupling to matter of heavy quarks is smaller than for 
light quarks as well as that the thermalization time for heavy quarks is longer than for light 
quarks. Due to the above features, one should observe a pattern of gradually increasing 
Raa when going from the mostly gluon-originated light-flavor hadrons {h^ and 7r°) to D 



to B mesons: i?^^ < i?^^ < i?^^ [2l|] (gluons lose more energy than quarks since gluons 
have a higher color charge). The enhancement above the unity of the heavy-to-light ratio 
^AA ~ ^aa/ -^AA probes the color charge dependence of the parton energy loss while the 
ratio R^A^ = Raa/ ^AA probes the mass dependence of the parton energy loss. 

Let us now consider what the experiment tells us about the open heavy flavor sup- 
pression and elliptic flow. As mentioned earlier, at RHIC, open heavy flavors can be studied 
through the measurements with the electrons and positrons originating from the semilep- 
tonic decays of D and B mesons. Fig. 10 shows the nuclear modification factor and elliptic 
flow coefficient for HF electrons as functions of pt, obtained in central collisions of gold 



nuclei at RHIC (closed circles) 22|- 25[. In contrast to the above expectations, the results 
for R^^ show a strong suppression of HF decay electrons at > 2 GeV/c, approaching 
at high Pt the level of suppression for 7r°. This evidences that produced medium is quite 
dense for heavy quarks to lose energy as efficiently as light quarks do. The measurement of 
elliptic flow gives rather large value for V2^. This means that HF electrons are involved in 
a collective motion being indicative of the collective flow of their parent particles as well. 

In Fig. [ini the results of some model considerations are also shown. The best description is 
provided by the model assuming the Brownian motion of heavy quarks within the framework 



of Langevin dynamics 26|. Let us consider the basic assumptions of this theory. Firstly, 



14 



DC 1.6^,1 (a) o-io%centrai Armesto et al. (I) 

1.4MIIIIIIIIII I I vanHeesetal. 

1.2: 



Au+Au @ vji;;;; = 200 Gev 

u. 0.2 J I I I I I I I I I I I I I I 

"i." : (b) 




7 8 9 

[GeV/c] 

FIG. 10: (Color online) Nuclear modification factor (upper panel, central Au-Au) and elliptic flow 
(lower panel, minimum-bias Au-Au) of non-photonic electrons at RHIC, compared to theory. The 
band corresponds to the Langevin simulations based on an expanding fireball with effective heavy 
quark resonance interactions 26 1 . 



the thermal heavy quark momentum ~ mT (in nonrelativistic approximation) is much 
larger than the typical momentum transfer ~ from a thermal medium to a heavy 
quark. Hence, the motion of a heavy quark in the QGP can be represented as the Brownian 
motion, which can be described by using the Langevin equation. Secondly, heavy quark 
loses its energy in elastic scattering processes with light partons. Besides, as evidenced from 
calculations of heavy and light meson correlators within LQCD, in QGP the D- and B- 
meson like resonant states exist up to T < 2Tc. Rescattering on these resonant states plays 
an important role in thermalizing heavy quarks. Thirdly, in order to get the spectrum of HF 
electrons, c- and 6-quarks are to be hadronized to D and B mesons via quark coalescence 
(at low pt) and fragmentation (at high p^). 

The analysis shows that resonance scattering decreases nuclear modification factor 
and increases azimuthal asymmetry . Heavy-light quark coalescence in subsequent 
hadronization significantly amplifies and increases -R^a ) especially in the pt — 2 GeV/c 
region. The contribution from B mesons to '^2^ estimated by providing full calcu- 
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FIG. 11: (Color online) Transverse momentum dependence of the relative contribution from B 
mesons to the non-photonic electron yields. The solid curve illustrates the fixed-order-plus-next- 
to-leading-log (FONLL) calculation 



lations with c+b quarks and with only c quarks. The result is that the B-meson contribution 
increases -R^a ^^"^ decreases V2^, and becomes important above — 3 GeV/c. The last 
point was actually confirmed in the recent PHENIX measurements of the bottom fraction of 
HF electrons, {b — )■ e)/[(6 — )■ e) + (c — )■ e)], obtained in p-p collisions at ^/s = 200 GeV and 
shown in Fig. 11 Thus, one can conclude, that the combined effect of coalescence of 
heavy quarks Q with light quarks q, and of the resonant heavy-quark interaction is essential 
in generating strong elliptic flow of up to 10%, together with strong suppression of 
heavy flavor electrons with about 0.5. 



IV. CHARMONIUM PRODUCTION AT FINITE TRANSVERSE MOMENTUM 



Now we can consider the production of J/ip mesons at finite transverse momentum pt 
and compare it with the corresponding dependence for open heavy flavor mesons. Fig . 12 



shows the pr-dependence of J/ip^s Raa in — 20% Cu+Cu collisions from PHENIX (29| and 



STAR 



30|, and — 60% Cu+Cu collisions from STAR 30|. The most striking peculiarity is 



that the nuclear modification factor for J /ip increases with px from the value of about 0.5 at 
low Pt to the value slightly exceeding unity at high > 5 GeV/c. These data indicate that 
there is no J /ip suppression at high pt in heavy-ion RHIC experiments - the result which 
could be considered as surprising taking into account the strong suppression of heavy flavor 
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FIG. 12: (Color online) J/V''s Raa vs. pt, obtained at RHIC experiments and compared with 
different model calculations. The best description of the data is provided by the two-component 
model [33[ (dotted line). 



electrons considered in the previous section. 

For comparison, Fig. 12 also shows the results of the model calculations for the open 
charm Raa- The solid line corresponds to the charm quark energy loss because of the elastic 
scattering and radiative parton processes, with assumed medium gluon density dNg/dy = 



254 for — 20% Cu+Cu collisions 3l|- The dash-dotted line shows the results of the 
model calculations for the D-meson energy loss with dNg/dy = 275 32|], where the D-meson 



suppression is caused by the coUisional dissociation in quark-gluon plasma. Both models, 
which correctly describe open heavy-flavor suppression in Au-Au collisions, predict charm 
meson suppression of a factor about 2 aX px > 5GeV/c, in contrast to the J/ip^s Raa- 
This comparison suggests that high-pj- J/tp production does not dominantly proceed via a 
channel carrying color, but rather the contribution of the color singlet channel is prevailing. 
The dotted line in Fig. [12] is the result for the J/ip's Raa{pt) in the two-component 



model [33[, which provides quite good qualitative description of the data. Let us consider 
the basic assumptions of this model in more detail. First, it assumes that the pt spectra of 
charmonia states (\& = J/tp,Xc,'ip') consist of two parts, direct and coalescence: 



dN^ dNq, 



PrdpT PrdpT 



_^ dN^ 
d^r PrdpT 



(5) 



coal 



The direct component is associated with hard production of charmonia in primordial 
N-N collisions, subject to suppression in the subsequent medium evolution due to gluon 
dissociation reactions in QGP and break-up by vr and p mesons in a hadron gas phase. The 
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FIG. 13: (Color online) Tevatron and STAR data for the ratio of J/^ from B-meson feed-down to 
inclusive J/'i/'- 

phase space distribution of different charmonia follows the Boltzmann transport equation, 
where the nuclear absorption is included in the initial conditions. Besides, the Cronin effect, 
consisting in multiple scattering of an initial parton on elementary constituents of the facing 
nucleus and leading to an increased pt in its final state (before fusion to charmonia), is 
also taken into account through initial conditions for charmonia momentum distributions. 
The model consistently treats the leakage effect, i.e., charmonia travelling outside the fire- 
ball boundary before freeze-out are not subject to dissociation. The leakage effect reduces 
suppression primarily for high-pj' charmonia. The soft component is associated with the 
coalescence of c, c quarks near the QCD phase boundary assuming an approximate thermal- 
ization up to c-quark momenta of pt ~ 2 — 2.5 GeV/c. The medium evolution is modelled 
by the isentropically expanding fireball with a cylindrical volume. 

Although in such a formulation the model is already quite complicated, two additional 
effects should be taken into account in order to reproduce an increasing trend of J/ip^s Raa 



at high Pt 



33|. These two additional effects are: 1. Finite formation time, required to 



build up the charmonium wave function from a " pre-hadronic" c-c pair - this effect leads to 
the reduction of the charmonia dissociation cross sections. 2. Feed-down from 5-mesons, 



Fig. [13] (from Ref. 33l|) shows the data from Tevatron 3^ and STAR on the B ^ J/iJj 
feed-down fraction in elementary p-p (p-p) collisions which is quite considerable at high pt- 
After including the formation time effects and B-meson feed-down (according to Tevatron 
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data) into the two- component model, it is able to reproduce the increasing behavior of 
J/ip^s Raa at high pt- Nevertheless, note that the recent STAR measurements on J/ip- 
hadron azimuthal correlations in p-p collisions at y/s = 200 GeV show that the J/ tp fraction 
from B-meson feed-down at high pt is not so significant, about 13% at > 5 GeV/c. This 
means that further studies are necessary, both theoretical and experimental, in order to 
provide a consistent picture of J ftp production at high p^. 
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